Effective Coulomb interactions within BEDT-TTF dimers 
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We calculate the effective Coulomb interactions between holes in dimers of the organic molecule 
BEDT-TTF in vacuo. We use density functional theory (DFT) to parameterise Hubbard models for 
j3 and k phase organic charge transfer salts. We focus on the intra-dimer Coulomb repulsion, , 
and the inter-monomer Coulomb repulsion, V&\ We find that U^ v) = 3.22 ± 0.09 eV and = 
2.71 ± 0.10 eV for 23 experimental geometries taken from a range of materials in both the /3 and k 
polymorphs. The quoted error is one standard deviation over the set of conformations studied. We 



conclude that Ujj" and Vm^ are essentially the same for an isolated dimer with the geometries present 
in all of the compounds studied. We highlight the disagreement between our parameterisation of 
the Hubbard model and previous results from both DFT and Hiickel methods and show that this is 



caused by the failure of an assumption made in previous calculations (that Um 3> Vi , where U, 
is the effective intra-monomer Coulomb repulsion). We discuss the implications of our calculations 
for theories of the BEDT-TTF salts based on the Hubbard model on the 2D anisotropic triangular 
lattice and explore the role of conformational disorder in these materials. 
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I. INTRODUCTION 



Layered organic charge transfer salts 
of the form (ET) 2 X, where ET is 
bis(ethylenedithio)tetrathiafulvalene or BEDT-TTF 
and X is a monovalent anion, exhibit a variety of 
unusual phenomena due to the strong electronic corre- 
lations present in these materials.— These phenomena 
include unconventional superconductivity^ with a small 
superfluid stiffness^ a Mott insulator^ a spin liquid, 4 
strongly correlated 5 ^ and unconventional 6 -^ metallic 
states, and a pseudogapJ* Experimentally, one can tune 
between these phases by varying the temperature and 
pressure (including both hydrostatic and 'chemical' 
pressure, i.e., varying the anion, X)^ 

DFT, as implemented with current approximate 
exchange-correlation functionals, does not capture sev- 
eral important aspects of the physics of strongly- 
correlated electronic systems. 9 For example, DFT band 
structure calculations of ET crystals , 10 ' 11 ' 12 produce a 
half-filled valence band and hence a metallic state. How- 
ever, these calculations do not recover the Mott insulat- 
ing state or the other strongly correlated effects that are 
observed experimentally. Therefore, efforts have focused 
on the application of effective low-energy Hamiltionians, 
such as Hubbard models J However, in molecular crys- 
tals, the effective parameters for such low energy Hamil- 
tonians may be calculated from studying the properties of 
single molecules or small molecular clusters, which may 
be accurately described by BFT^^^Ll^. 

ET salts occur in a variety of crystal packing struc- 
tures. In the /? and k polymorphs the ET molecules form 
dimers. Intradimer dynamics are often integrated out of 
effective low energy models of /?-(ET) 2 X and k-(ET) 2 X. 
In these charge transfer salts, each dimer donates one 
electron to the anion layers to form a half-filled system. 
Both Hiickel^ and DFTi * 1 !^ calculations have found 



that the dimers form an anisotropic triangular lattice in 
which each lattice site is a dimer. However, there is a 
strong effective Coulomb repulsion, Ud, between two elec- 
trons on the same dimer, which must be included in the 
effective low energy description, 1 "^ 1 - The electronic inter- 
actions within an ET dimer are stronger than those be- 
tween an ET molecule and its next-nearest neighbours 
on the crystal lattice. Therefore, these materials have 
been widely studied on the basis of Hubbard models. 1 
In order to explain the observed physics these theories 
assume that both chemical and hydrostatic pressure re- 
duce Ud/W, where W is the bandwidth. Therefore, an 
important task for the field is to understand how this 
ratio varies with chemical and hydrostatic pressure. 

Previously, the on-site Coulomb repulsion term in 
the Hubbard model, Ud, has been estimated from both 
Hiickel^^^^^^^^ and density functional^ 
calculations under the assumptions that the intra- 
monomer Coulomb repulsion U m — > oo, and the inter- 
monomer Coulomb repulsion V m — > 0. We will show 
below that this assumption is incorrect and leads to a 
systematic underestimate of Ud- 

Disorder plays a number of important roles in or- 
ganic superconductors. 31 Increasing the degree of disor- 
der leads to a suppression of the superconducting critical 
temperature, T c , which is correlated with a rise in the 
residual resistivity. 31 Further disorder can cause a viola- 
tion of Matthiessen's rule via impurity assisted tunnelling 
in the interlayer direction— 

In K-(ET) 2 Cu[N(CN) 2 ]Br the degree of disorder can 
be increased by increasing the rate at which the sam- 
ple is cooled ) 33 ' 34 ' 35 ' 36 ' 37 which leads to a suppression 
of T c by ~ 1 K. Further, increasing the cooling rate 
can drive the system towards the insulating side of the 
metal-insulator transition^ Two hypothesis have been 
proposed for the source of this disorder: terminal ethy- 
lene group disorde r 33 ' 34 ! 35 ' 36 and disorder in the anion 
layer^ Therefore it is important to estimate the scat- 



tering rate caused by terminal ethylene disorder in this 
material. 

Even more dramatic effects are observed in /3-(ET)2l3- 
Variations of the pressure as the sample is cooled can 
change the ambient pressure T c from 1 K (for samples 
cooled at ambient pressure; known as the (3l phase) to 7 
K (for samples cooled at P > 1 kbar once the pressure 
is released; known as the (3h phase) £2- In this material 
clear differences in the terminal ethylene groups in the 0h 
and Pl phases are observed via x-ray scattering^ Thus 
it has been argued that the terminal ethylene disorder is 
responsible for the differences in the critical temperatures 
between the (3h and (3l phases4i 

Therefore, we also present calculations of the effective 
site energy for holes, £,j> f° r the f3- and k- phase salts. 
This allows us to study the effects of impurity scattering 
caused by conformational disorder in the terminal ethy- 
lene groups of the ET molecule. 

In this paper we present DFT calculations for ET 
dimers in vacuum. In Sec. II we describe the compu- 
tational method by which we calculate these energies. 
In Sec. Ill we discuss the problem for the isolated dimer 
and review the parameterisation of the two-site extended 
Hubbard model from the total energies of the relevant 
(ET) 2 charge states. In Sec. IV we report and discuss 

(v) (v) 

the resulting values of U d , Vm and £d- In Sec. V we 
draw our conclusions. 



II. COMPUTATIONAL METHODS 

We used DFT to calculate the total energies of ET 
dimers in various conformations and charge states. We 
used the SIESTA 42 implementation of DFT, with the 
PBE exchange-correlation functional, 43 a triple-^ plus 
single polarisation (TZP) basis set (except where we ex- 
plicitly specify otherwise) and basis functions consisting 
of Sankey type numerical atomic orbitals^ 4 - The orbital 
functions were confined to a radius r c from their cen- 
tres, which slightly increases the energy of the orbital. 
The specified maximum allowed increase in energy due 
to this cutoff was 2 mRy. The convergence of the in- 
tegration mesh was determined by specifying an effec- 
tive plane-wave cutoff energy of 250 Ry. The initial 
spin moments on each atom were arranged antiferromag- 
netically wherever possible. We used pseudopotentials 
constructed according to the improved Troullier-Martins 
(TM2) method^ 

Nuclear positions for C and S atoms were obtained 
from x-ray crystallographyi 22 i 23 i 46 i 47 i 48 i 49 i 50 ' 51 H atoms, 
which are not observed in x-ray scattering experiments, 
were relaxed by the conjugate-gradient method. To- 
tal DFT energy differences between the relevant charge 
states [E(l) - E(0)} and [E(2) - E(0)] were equated 
with the corresponding analytical expressions in Eq. ([5]) 
to determine the Hubbard parameters. We focus on these 
'experimental' geometries rather than performing a full 
relaxation for a number of reasons. Firstly, there are 



small differences in the reported geometries for different 
ET salts, and one would like to understand the effect 
of these. Secondly, the experiments effectively 'integrate 
over' all of the relevant charge states and therefore pro- 
vide an 'average' conformation. Thirdly, the experiments 
naturally include the effects on the molecular conforma- 
tion due to the crystalline environment, which are absent 
from in vacuo calculations. 



III. THE TWO SITE EXTENDED HUBBARD 
MODEL 

In calculations of the effective Coulomb in- 
teraction (the Hubbard U) in molecular solids 
it is important to recognise that their are two 
contributions^^ 1 ^! 1 -^ 1 ^ That is, the effective 
Coulomb interaction on a ET dimer may be written as 



where U d is the value of Ud for the dimer cluster in vac- 

(p) 

uum, and U d is the reduction in Ud from the polarisable 

(%>) 

crystalline environment. Calculating U d for ET salts is 
a highly non-trivial problem due to the large size of the 
ET molecule relative to the intermolecular spacing. Be- 

(v) 

low we present results of DFT calculations for U d of 
dimers in the conformations found in a wide range of k 
and j3 phase ET salts. Similar results hold for U m , the ef- 
fective Coulomb repulsion between two holes on the same 
monomer and V m the effective Coulomb interaction be- 
tween two holes on neighbouring monomers. Below we 
will primarily discuss the vacuum contributions to these 

terms, Um and Vm . 

The effective vacuum intradimer Coulomb energy, 

Jjfj" is given by (see, e.g., Ref. [H) 

U^=E(0) + E(2)-2E(1), (2) 

where E(q) is the ground state energy of the dimer in vac- 
uum containing q holes, i.e., with charge +q. Similarly, 
the effective site energy for holes is given by 

£ d = E(Q) - E(l). (3) 

Below we calculate E(q) via density functional methods. 

It is also interesting to consider intradimer dynamics, 
which can be described via a two site extended Hubbard 
model r 

H = ^ ^rnSia- - t ^ (h\ a h 2<T + h.c)j 

+ ^ Umifiiihii + V m hin 2 (4) 

i 

where h^J annihilates (creates) a hole on site (monomer) 
i in spin state cr, £ m i is the site energy for holes on site 
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FIG. 1: The HOMO of (ET)% + (top) and charge neutral 
(ET)2 (bottom), with nuclear positions from the crystal /3- 
(ET) 2 I 3 . The HOMO of (ET)^ + is the dimer bonding orbital 
and the HOMO of (ET) 2 is the antibonding orbital of the two 
ET HOMOs (cf. Fig. 0. The essential difference between the 
two lies in the relative phase of the orbital function on each 
molecule. The bonding orbital connects the ET molecules at 
the S- ■ ■ S contacts (cf. Fig. In the antibonding orbital, 
there are nodes between the S- • • S contacts. 

i, flic is the number operator for spin a holes on site i, 
hi = J2a ^io-j t is the intradimer hopping integral, U m i is 
the effective on-site (monomer) Coulomb repulsion, and 
V m is the intersite Coulomb repulsion. 

The lowest energy eigenvalues of Hamiltonian Q for 
each charge state are 

E(0) = 0, (5a) 

E(l) = ? m -iv/4*2 + (A£ m )2, (5b) 

and 

E(2) = 2l m + ~ (2U m + V m - 2 A cos 9) (5c) 

where £ m = i(£ mX + £ m2 ), A = \2t 2 + {AU m ) 2 + (U ml - 
V m )(U m 2 - V m ) + 3(A^ m ) 2 , cos39 = (U m - 2V m )[18t 2 - 
(2U ml - U m2 - V m )(U ml - 2U m2 + V m ) - 9(A^ m ) 2 ]/2A 3 . 

^Cm = Cml — S,m2, U m — ^(Uml + Um2), and AU m = 
U m l — U m 2- 



FIG. 2: The HOMO of (ET)^ + (top) and charge neutral 
(ET) 2 (bottom), with nuclear positions from the crystal k- 
(ET)2Cu2(CN)3. The similarity of the nuclear structures and 
orbitals between this conformation and the /3 conformation 
in Fig. [T] highlight the dimer as a common structural unit 
within two different packing motifs. 




FIG. 3: The HOMO of a charge neutral ET monomer, with 
nuclear positions from the crystal /3-(ET)2l3. This is the or- 
bital from each molecule that contributes to the HOMO of 
the (ET) 2 and (ET)^+ dimers. 



We have previously calculated £ m and U m from the 
one site Hubbard model for an ET monomer for the ex- 
perimental observed conformations in all the materials 
discussed below therefore one may solve Eqs. ([5]) for t 

(v) 

and V m taking £ m and U m from the monomer calcula- 
tions. The case of two holes with different site energies 
and on-site Coulomb repulsion may be solved by a general 
method for diagonalising cubic matrix eigensystems^ In 
cases where the two molecules in a dimer have the same 
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geometry (e.g., by symmetry), £ ml = £ m2 = £ m and 
U m i — U m 2 = U m and the eigenvalues simplify to 

E(l) = U-t (6a) 

E(2) = 2U + 1 (ll m + V m - v/l6t 2 + (U m - F m ) 2 ) 

(6b) 

in which case the solution is straightforward. 

In the limit U m = V m = the two site Hubbard model 
has two solutions: the bonding state \4>ba) = |01ct} + |</ , 2ct) 
and the antibonding state \4>ba) = l^ier) — \4>2a)i where 
\4>ia) = hicr\0) is a single electron state centred on the i th 
monomer and |0) is the (particle) vacuum state. 

In Figs. CD and [2] we plot the HOMOs of ET 
dimcr for the conformations found in /3-(ET)2l3 and k- 
(ET) 2 Cu2(CN) 3 respectively, in both the charge neutral 
and the 2+ states. It can be seen that these dimer or- 
bitals are the antibonding and bonding hybrids of the 
ET monomer HOMO (shown in Fig. [3]), respectively. 
The most important difference between the orbital ge- 
ometries lies in the S- • • S intermolecular contacts, which 
contain nodes in the antibonding orbital, but are con- 
nected in the bonding orbital. Thus the DFT picture of 
the (ET)2 system is remarkably similar to the molecular 
orbital description of a diatomic molecule^ but with the 
'covalent bond' between the two monomers rather than 
between two atoms. 




FIG. 4: Variation of Hubbard parameters found from DFT 
calculations with basis set, which improve from left to 
right. The test conformation is taken from the crystal k- 
(ET) 2 Cu(NCS) 2 .— All of the quantities except t are well- 
converged at TZP, the basis set chosen for all subsequent cal- 
culations. Indeed, £ m is the only other quantity that changes 
significantly across the range of basis sets. However, t is a 
relatively small quantity, on the order of its own variations 
with respect to basis set size. Hence we conclude that solving 
the Hamiltonian @ is not an accurate method for finding t. 



IV. CALCULATION OF THE HUBBARD 
MODEL PARAMETERS 

A. Basis set convergence 

We tested the basis set convergence of the DFT 
calculations using the conformation observed in re- 
(ET) 2 Cu(NCS) 2 as the test case, with single-C (SZ), 
single-^ plus polarisation (SZP), double-^ (DZ), double- 
£ plus polarisation (DZP) and TZP basis sets. We also 
calculated the monomer parameters, U m and £ m , in each 
basis, using the method we previously applied to the ET 
monomers.— The Hubbard model parameters in each ba- 
sis set are reported in Fig. 0J The values of all parameters 
are well- converged in the TZP basis, except t. t is an or- 
der of magnitude smaller than the other parameters, and 
on the order of both the variation of the other param- 
eters among the basis sets tested and the uncertainty 
associated with the calculation method. This suggests 
that extracting t from band structure calculation a 11 ' 12 is 
a more accurate and reliable method of estimating the 
hopping integrals in these systems. 
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FIG. 5: The effective intra-dimer Coulomb repulsion, , for 
various ET dimers. The a;-axis separates the data by source 
crystal polymorph (/3 or «), and by the terminal ethylene 
group conformation of each ET molecule in the dimer. Uj V ' 
does not change significantly across the different ET crystals 
examined. For /3-(ET) 2 X crystals, U { d v) = 3.19 ±0.07 eV. For 



B. Variation of the intra-dimer Coulomb repulsion k-(ET) 2 X crystals, U 
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3.23 ± 0.09 eV. The difference in 



Now we consider variation of the Hubbard model pa- 
rameters across the conformations found in different ma- 
terials, beginning with In Fig. Owe show the values 



between the two crystal polymorphs is ~ 1%. There- 
fore, there is no significant dependence of on the dimer 
geometry associated with different crystal polymorphs. 
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FIG. 6: ET molecules within the crystals studied occur in 
two conformations, denoted eclipsed and staggered. The dif- 
ference between them lies in the relative orientation of the 
terminal ethylene groups. 



of l/jj"' 1 for the conformations observed experimentally in 
a variety of ET crystals. Of particular note are the three 
data points corresponding to different possible confor- 
mations of /3-(ET)2l3. In the ET molecule the terminal 
ethylene groups may take two relative orientations known 
as the staggered and eclipsed conformations (cf. Fig. [6]). 

(v) 

U d is smallest when both ET molecules are in the stag- 

(v) 

gered conformation. Conversely, the largest U d for this 
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FIG. 7: Dimer hole site energy, £d, for various ET dimers. 
For /3-(ET) 2 X crystals, the mean value is £ d = 4.45 ± 0.10 
eV. For k-(ET)2X crystals, the mean value is £d = 4.46 ±0.14 
eV. The mean value for the whole data set is £d = 4.45 ± 0.13 
eV. The effect of conformation on ^ is significantly larger for 
/3-(ET)2l3 (~ 10%) than it is for the other parameters. The 
variations of £d with dimer geometry associated with crys- 
tal polymorph and anion are ~ 3%, similar to the relative 



crystal occurs when both ET molecules are eclipsed, with variations of and Vm' across the whole data set 
intermediate U d values for the case with one staggered 
and one eclipsed ET molecule. This trend is repeated in 
the K-phase crystals, where two data sets (corresponding 
to different temperatures at which the nuclear positions 
were determined) for k-(ET)2Cu[N(CN)2]I provide data 
for both conformations. 

The mean value of U d for the p phase crystals is 3.19 

(v) 

± 0.07 eV, and the mean U\ for the k phase crystals is 
3.23 ± 0.09 eV. The quoted error is one standard devia- 
tion over the full set of conformations studied. The dif- 
ference between the two values of U d is ~ 1%, and well 

' (v) 

within the error ranges. This suggests that U d takes 
the same value, 3.22 ± 0.09 eV, in all (3 and k phase ET 
salts. This result is significantly larger than the value of 
U d obtained from Hiickel calculations (~0.5-2 eV), as 
we will discuss below. 



r(v) 



C. Variations in site energy and the role of disorder 

As reviewed in the introduction, a number of exper- 
iments have shown that disorder has strong effects on 
both the normal state and superconducting properties 
of organic charge transfer salts There 
has been relatively little work on the effect of the ran- 
dom U Hubbard model. Conclusions drawn from studies 
in one dimensio n 53 i 54 cannot be straightforwardly gen- 
eralised to higher dimensions. Mutou» used dynami- 
cal mean field theory to study the metallic phase of the 



random U Hubbard model. However, he did not con- 
sider the effect a random U on either superconductivity 
or the Mott transition, which are the primary concerns in 
the organic charge transfer salts. However, Mutou con- 
cluded that for small impurity concentrations Kondo-likc 
effects mean that the random U Hubbard model is signifi- 
cantly different from the virtual crystal approximation to 
the random U Hubbard model, which describes the sys- 
tem in terms of an average U. The only study^ we are 
aware of that discusses superconductivity in the random 
U Hubbard model treats the negative U model, which is 
not realistic for the organic charge transfer salts. Litak 
and Gyorffy^i studied a model where some sites have 
U = and others have a negative U. They find that 
superconductivity is suppressed above at certain critical 
concentration of U = sites. Therefore, it is not clear 
what implications our finding of small changes in Uj" 
and hence Ud has for the physics of the organic charge 
transfer salts. However, it is interesting to ask what role 
this plays in the observed role of disorder in suppressing 
superconductivity^ and driving the system towards the 
Mott transition^ 

To understand the role of conformational disorder in 
terms of an effective Hamiltonian built up from ET 
dimers one must also understand the effect of confor- 
mational disorder of the effective dimer site energy (for 
holes), £<j. This is straightforwardly found from the DFT 
calculations described above via Eq. © and the results 
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FIG. 8: Intradimer for various ET dimers. For /3-(ET) 2 X 



crystals, = 2.69 ± 0.13 eV and for k-(ET) 2 X crystals, 



V, 
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2.72 ± 0.09 eV. The mean value is V, 



(v) 



2.71 ±0.10 



eV. The difference in Vrk between the crystal polymorphs 
is ~ 2%. Therefore, Vm \ like Ufj°\ does not significantly 
depend on the geometry associated with crystal polymorph. 

(v) 

The effect of ET conformation on the value of Vm in the 
crystals /3-(ET) a Ia and k-(ET) 2 Cu[N(CN) 2 ]I is also similar 
to the effect on Uj" . Vm is lowest when the ET dimer has 
the staggered-staggered conformation, and rises when either 
or both ET molecules are eclipsed. 



are reported in Fig. [71 The effective scattering rate due 
to conformational disorder is given by 



= J2N i 7rD(E F )\A i U\ 



(7) 



where i labels the type of impurity (both staggered or 
mixed; the ground state conformation is both eclipsed), 
Ni is number of impurities of type i, D(Ep) is the density 
of states at the Fermi level, and Aj£d is the difference 
between ^ for i type impurities and ^ of eclipsed dimers. 

In quasi-2D systems, D(E F ) is simply related to the 
cyclotron electron mass^i by the relation 



D{E F ) 



m c 
2nh 2 



(8) 



and in the presence of interactions Luttinger's theorem 58 
for a Fermi liquid produces 



D{E F ) 



2ttH 2 



(9) 



where m* is the effective mass. From quantum os- 
cillation measurements, Wosnitza et aL 5 = found that 
m*/m e — 4.2 in /3-(ET) 2 I 3 , where m e is the electron 
rest mass. From Shubnikov-de Haas measurements in k- 
(ET) 2 Cu[N(CN) 2 ]Br Caulfield et dM found that m*/m e 



= 6.4. The scattering rate r can be found from mea- 
surement of the interlayer residual resistivity, po , by the 
relation 61 



Po 



2e 2 m*ct 2 L T 



(10) 



where c is the interlayer lattice constant taken from 
the relevant x-ray scattering measurements* 6 -^ and t± 
is the interlayer hopping integral, which has previously 
been estimated from experimental data for both k- 
(ET) 2 Cu[N(CN) 2 ]Br (Ref. M) and /3-(ET) 2 I 3 (Ref. HH). 
Using these parameters we calculated the scattering rate 
in both the (3l and (3h phases of /3-(ET) 2 l3 from the 
low temperature values of po reported by Ginodman et 
al.£& The scattering rate due to conformational impu- 
rities, t" 1 is then T7 1 = r 



H 



-£ , where r L (t h ) is 
the quasiparticle lifetime in the (3l {J3h) phase. Given 
our calculated values of Ai£d an ~ 8% concentration 
of staggered impurities would be required to cause this 
scattering rate. From a similar calculation comparing 
the residual resistivity measured in a single sample of 
K-(ET) 2 Cu[N(CN) 2 ]Br cooled at different rates we find 
that a ~ 2% concentration of staggered impurities would 
be sufficient to explain the rise increase in residual re- 
sistivity observed in the experiment utilising the fastest 
cooling over that performed with the slowest cooling rate. 
X-ray scattering experiments^ 8 - find that 3 ± 3% of the 
ET molecules are in the staggered conformation at 9 K, 
which is entirely consistent with our result. However, 
Wolter et aZ.'s 38 argument that this impurity concentra- 
tion is too small to cause the observed effects of disorder 
in not sustained by the above calculations. Rather we 
find that all of the suppression in T c is entirely consis- 
tent with this degree of disorder. 



D. Variations in inter- molecular Coulomb repulsion 

In Fig. [8] we show the distribution of the calculated 
values of vJk ] ■ The mean value of Vm for the [3 phase 
crystals is 2.69 ± 0.13 eV, while the mean value for the 
K phase crystals is 2.72 ± 0.09 cV. Again, the difference 
between the values is small (~ 2%) compared to the dis- 
tribution for each polymorph. Therefore, Vm^ is essen- 
tially the same across all of the conformations studied, 
with a mean value of 2.71 ± 0.10 eV. 

(v) 

Previous calculations of U) based on both the Huckel 
methocl^^^^^S^^^ and DFT— ^ have assumed 
that Um^ oo and Vm — 0. Substituting these condi- 
tions into Eqs. ^ and © yields = 2t. Literature 

values of based on this approximation are presented 
in Table U for comparison with our DFT results. It can 
be seen that this assumption yields values of that 
are significantly smaller than those we have calculated 
above (cf. Fig. 0. However, we have previously founcU^ 

(v) 

that U m = 4.2 ±0.1 eV. Comparing this with the above 
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Crystal 


Method XJ[ 


i V) (eV) 


/3-(ET) 2 I 3 


Hiickel 


0.49 


/3-(ET) 2 IBr 2 


Hiicke! 25 


0.98 


/3-(ET) 2 ICl 2 


Hiickel 


1.04 


/3-(ET) 2 I 3 


Hiickel 2 ^ 


0.49 


/3-(ET) 2 CH(S0 2 CF 3 ) 2 


Hiickel 0. 


88-0.90 


/3-(ET) 2 [OsNOCl 5 ] 


Hiickel 


2.10 


k-(ET) 2 Cu[N(CN) 2 ]C1 


DFT 12 


0.4 


K-(ET) 2 Cu[N(CN) 2 ]Br 


Hiickel 


0.45 


k-(ET) 2 Cu(NCS) 2 


Hiicke! 29 


0.48 


k-(ET) 2 Cu(NCS) 2 


Hiicke! 23 


0.14 


k-(ET) 2 Cu(NCS) 2 


Hiicke! 30 


0.46 


k-(ET) 2 Cu(NCS) 2 


DFT— 


0.83 


k-(ET) 2 Cu(NCS) 2 


DFT 12 


0.4 


K-(ET) 2 Cu[N(CN) 2 ]Br 


Hiicke! 3 " 


0.49 


k-(ET) 2 Cu 2 (CN) 3 


Hucke! 30 


0.45 


k-(ET) 2 Cu 2 (CN) 3 


DFT— 


0.85 




DPT 12 


n d 


«-(ET) 2 I 3 


Hiickel 


0.49 


«-(ET) 2 I 3 


Hiicke^ 


0.22 


I: Previous estimates 


of U { d v) for 


various 



K-phase ET salts. These values were obtained from both 
Hiickel and density functional methods under the assump- 
tions Um oo and Vm = 0, which yields U d = It. These 
estimates substantially underestimate the actual value of l/jj"^ 

(see Fig. [5]) as ~ Vm • The two site extended Hubbard 
model produces values of t on the same order of magnitude 
as these Hiickel calculations. One should also note the wide 
scatter between the different Hiickel calculations, even be- 
tween different studies of the same material. 

results we see that Um /Vm ~ 1-5, in contradiction to 
the assumption that Um ^> Vm • Hence Um 3> \2t\. 

If we instead make the assumption Um — Vm 3> 
then Eqs. © and © give 

u^^liu^ + vW). (11) 

Substituting in the mean values of Um) and Vm^ gives 
U d = 3.41 eV. This result is close (within 6%) to our 

(v) 

calculated value of U d . Therefore, this is a reasonable 
approximation for the ET salts. Further, this shows that 

(v) 

the result that U d does not vary significantly because of 
changes in conformation between different salts or poly- 
morphs is a consequence of the fact that neither Um or 
Vm vary significantly because of changes in conforma- 
tion between different salts or polymorphs. 

V. CONCLUSIONS 

The effective Coulomb repulsion terms in the Hubbard 
model are essentially the same for all of the ET con- 
formations studied. We found that U ( d v} = 3.22 ± 0.09 

eV, V}n ] = 2.71 ± 0.10 eV. The value of U d v) is sig- 
nificantly larger than previous estimates based on the 



extended Hiickel formalism or DFT under the assump- 
tions Um — * 00 and Vm'^ = 0. This can be under- 

(v) (v) 

stood because we have shown that Um ~ V m and hence 

The lack of variation of U d between the two poly- 
morphs and when the anion is changed is interesting in 
the context of theories of these organic charge transfer 
salts based on the Hubbard model. These theories require 
Ud/W to vary significantly as the anion is changed (chem- 
ical pressure) and under hydrostatic pressure. Therefore 

(r>) 

our results show that either UY or W must vary signif- 
icantly under chemical and hydrostatic pressure, or else 
these theories do not provide a correct description of the 
P and k phase organic charge transfer salts. This is par- 
ticularly interesting as fast cooling has been shown to 
drive K-(ET) 2 Cu[N(CN) 2 ]Br to the insulating side of the 
metal-insulator transition.— 

We have also studied the effects of conformational dis- 
order on these parameters, which is found to be quite 
small, consistent with the often subtle effects of con- 
formational disorder observed in these materials. The 
largest changes are found in the geometries taken from 
/3-(ET) 2 l3, which shows the strongest effects of confor- 
mational disorder. It is also interesting that we found a 

(v) 

systematic variation in U d is caused by conformational 
disorder. As there has been relatively little work on the 
random U Hubbard model it is difficult to speculate what 
effects this has on the low temperature physics of the or- 
ganic charge transfer salts at present. 

Given that DFT band structure parameterisations of 
the interdimer hopping integrals have recently been re- 
ported for several organic charge transfer salts ; 11 ^ 2 the 
outstanding challenge for the parameterisation of the 
Hubbard model in these systems is the accurate calcula- 
tion of l&\ The bandwidth in both the (3 (Ref. |T(J and 
k (Ref. [l"llfl2l ) phase salts is around 0.4-0.6 eV. Therefore, 
our finding that Uj"^ is significantly larger than has been 
realised previously shows that U d must be significant 
as if Ud — U^ then all of these materials would be well 
into the Mott insulating regime. Thus U d p ^ must signifi- 
cantly reduce Ud in order for the, observed, rich phase di- 
agram to be realised. This is consistent with comparisons 
of DMFT calculations to optical conductivity measure- 
ments on «-(ET) 2 Cu[N(CN) 2 ]Br a; Cli_ !1 ., which suggest 
that Ud — 0.3 eVi£2 Further, U^ may be quite sensitive 
to the crystal lattice and therefore may be important for 
understanding the strong dependence of these materials 
on chemical and hydrostatic pressure. 
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